Accelerated Computation of Free Energy Profile at ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semi-Empirical Reference Potential. I. Weighted Thermodynamics Perturbation

Free energy profile (FE Profile) is an essential quantity for the estimation of reaction rate and the validation of reaction mechanism. For chemical reactions in condensed phase or enzymatic reactions, the computation of FE profile at ab initio (ai) quantum mechanical/molecular mechanics (QM/MM) level is still far too expensive. Semiempirical (SE) method can be hundreds or thousands of times faster than the ai methods. However, the accuracy of SE methods is often unsatisfactory, due to the approximations that have been adopted in these methods. In this work, we proposed a new method termed MBAR+wTP, in which the ai QM/MM free energy profile is computed by a weighted thermodynamic perturbation (TP) correction to the SE profile generated by the multistate Bennett acceptance ratio (MBAR) analysis of the trajectories from umbrella samplings (US). The weight factors used in the TP calculations are a byproduct of the MBAR analysis in the post-processing of the US trajectories, which are often discarded after the free energy calculations. The results show that this approach can enhance the efficiency of ai FE profile calculations by several orders of magnitude

Representation of the QM Subsystem for Long-Range Electrostatic Interaction in Non-Periodic Ab Initio QM/MM Calculations

This paper is published as part of a thematic issue of Molecules on “Combined Quantum Mechanical and Molecular Mechanical Methods and Simulations”. http://www.mdpi.com/journal/molecules/special_issues/QMIn QM/MM calculations, it is essential to handle electrostatic interactions between the QM and MM subsystems accurately and efficiently. To achieve maximal efficiency, it is convenient to adopt a hybrid scheme, where the QM electron density is used explicitly in the evaluation of short-range QM/MM electrostatic interactions, while a multipolar representation for the QM electron density is employed to account for the long-range QM/MM electrostatic interactions. In order to avoid energy discontinuity at the cutoffs, which separate the short- and long-range QM/MM electrostatic interactions, a switching function should be utilized to ensure a smooth potential energy surface. In this study, we benchmarked the accuracy of such hybrid embedding schemes for QM/MM electrostatic interactions using different multipolar representations, switching functions and cutoff distances. For test systems (neutral and anionic oxyluciferin in MM (aqueous and enzyme) environments), the best accuracy was acquired with a combination of QM electrostatic potential (ESP) charges and dipoles and two switching functions (long-range electrostatic corrections (LREC) and Switch) in the treatment of long-range QM/MM electrostatics. It allowed us to apply a 10Å distance cutoff and still obtain QM/MM electrostatics/polarization energies within 0.1 kcal/mol and time-dependent density functional theory (TDDFT)/MM vertical excitation energies within 10−3 eV from theoretical reference values.This research was funded by the U.S. Department of Energy Office of Science (DE-SC0011297), EPSRC (EP/R013012/1, EP/L027151/1 and EP/N020669), and ERC (Project 757850 BioNet).Ye

DeepPredict : A zone preference prediction system for online lodging platforms

Publisher Copyright: © The author(s) 2021.Online lodging platforms have become more and more popular around the world. To make a booking in these platforms, a user usually needs to select a city first, then browses among all the prospective options. To improve the user experience, understanding the zone preferences of a user's booking behavior will be helpful. In this work, we aim to predict the zone preferences of users when booking accommodations for the next travel. We have two main challenges: (1) The previous works about next information of Points Of Interest (Pals) recommendation are mainly focused on users' historical records in the same city, while in practice, the historical records of a user in the same city would be very sparse. (2) Since each city has its own specific geographical entities, it is hard to extract the structured geographical features of accommodation in different cities. Towards the difficulties, we propose DeepPredict, a zone preference prediction system. To tackle the first challenge, DeepPredict involves users' historical records in all the cities and uses a deep learning based method to process them. For the second challenge, DeepPredict uses HERE places API to get the information of pals nearby, and processes the information with a unified way to get it. Also, the description of each accommodation might include some useful information, thus we use Sent2Vec, a sentence embedding algorithm, to get the embedding of accommodation description. Using a real-world dataset collected from Airbnb, DeepPredict can predict the zone preferences of users' next bookings with a remarkable performance. DeepPredict outperforms the state-of-the-art algorithms by 60% in macro Fl-score.Peer reviewe

Machine-Learning-Assisted Free Energy Simulation of Solution-Phase and Enzyme Reactions

Despite recent advances in the development of machine learning potentials (MLPs) for biomolecular simulations, there has been limited effort on developing stable and accurate MLPs for enzymatic reactions. Here we report a protocol for performing machine-learning-assisted free energy simulation of solution-phase and enzyme reactions at the ab initio quantum-mechanical/molecular-mechanical (ai-QM/MM) level of accuracy. Within our protocol, the MLP is built to reproduce the ai-QM/MM energy and forces on both QM (reactive) and MM (solvent/enzyme) atoms. As an alternative strategy, a delta machine learning potential (ΔMLP) is trained to reproduce the differences between the ai-QM/MM and semiempirical (se) QM/MM energies and forces. To account for the effect of the condensed-phase environment in both MLP and ΔMLP, the DeePMD representation of a molecular system is extended to incorporate the external electrostatic potential and field on each QM atom. Using the Menshutkin and chorismate mutase reactions as examples, we show that the developed MLP and ΔMLP reproduce the ai-QM/MM energy and forces with errors that on average are less than 1.0 kcal/mol and 1.0 kcal mol–1 Å–1, respectively, for representative configurations along the reaction pathway. For both reactions, MLP/ΔMLP-based simulations yielded free energy profiles that differed by less than 1.0 kcal/mol from the reference ai-QM/MM results at only a fraction of the computational cost

Pre‐symptomatic transmission of novel coronavirus in community settings

We used contact tracing to document how COVID‐19 was transmitted across 5 generations involving 10 cases, starting with an individual who became ill on January 27. We calculated the incubation period of the cases as the interval between infection and development of symptoms. The median incubation period was 6.0 days (interquartile range, 3.5‐9.5 days). The last two generations were infected in public places, 3 and 4 days prior to the onset of illness in their infectors. Both had certain underlying conditions and comorbidity. Further identification of how individuals transmit prior to being symptomatic will have important consequences.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/163478/2/irv12773.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/163478/1/irv12773_am.pd

TNFRSF10C methylation is a new epigenetic biomarker for colorectal cancer

Background Abnormal methylation of TNFRSF10C was found to be associated with different types of cancers, excluding colorectal cancer (CRC). In this paper, the performance of TNFRSF10C methylation in CRC was studied in two stages. Method The discovery stage was involved with 38 pairs of CRC tumor and paired adjacent non-tumor tissues, and 69 pairs of CRC tumor and paired adjacent non-tumor tissues were used for the validation stage. Quantitative methylation specific PCR (qMSP) method and percentage of methylated reference (PMR) were used to test and represent the methylation level of TNFRSF10C, respectively. A dual-luciferase reporter gene experiment was conducted to evaluate the promoter activity of TNFRSF10C fragment. Results A significant association of TNFRSF10C promoter hypermethylation with CRC was found and validated (discovery stage: 24.67 ± 7.52 vs. 3.36 ± 0.89; P = 0.003; validation stage: 31.21 ± 12.48 vs. 4.52 ± 1.47; P = 0.0005). Subsequent analyses of TCGA data among 46 pairs of CRC samples further confirmed our findings (cg23965061: P = 4E − 6; cg14015044: P = 1E − 7). Dual-luciferase reporter gene assay revealed that TNFRSF10C fragment was able to significantly promote gene expression (Fold change = 2.375, P = 0.013). Our data confirmed that TNFRSF10C promoter hypermethylation can predict shorter overall survival of CRC patients (P = 0.032). Additionally, bioinformatics analyses indicated that TNFRSF10C hypermethylation was significantly associated with lower TNFRSF10C expression. Conclusion Our work suggested that TNFRSF10C hypermethylation was significantly associated with the risk of CRC